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CHANNEL ESTIMATION FOR TIME DIVISION 
DUPLEX COMMUNICATION SYSTEMS 

This application claims priority to U.S. Provisional Patent Application No. 
60/175,167, filed on January 7, 2000. 

BACKGROUND 

The invention generally relates to wireless communication systems. In particular, the 
invention relates to channel estimation in a wireless communication system. 

Figure 1 is an illustration of a wireless communication system 10. The 
communication system 10 has base stations to 12 5 which communicate with user 
equipments (UEs) 14j to 14 3 . Each base station 12j has an associated operational area where 
it communicates with UEs I4 l to 14 3 in its operational area. 

In some communication systems, such as code division multiple access (CDMA) and 
time division duplex using code division multiple access (TDD/CDMA), multiple 
communications are sent over the same frequency spectrum. These communications are 
typically differentiated by their chip code sequences. To more efficiently use the frequency 
spectrum, TDD/CDMA communication systems use repeating frames divided into time slots 
for communication. A communication sent in such a system will have one or multiple 
associated chip codes and time slots assigned to it based on the communication's bandwidth. 

Since multiple communications may be sent in the same frequency spectrum and at 
the same time, a receiver in such a system must distinguish between the multiple 
communications. One approach to detecting such signals is single user detection. In single 
user detection, a receiver detects only the communications from a desired transmitter using 
a code associated with the desired transmitter, and treats signals of other transmitters as 
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interference. Another approach is referred to as joint detection. In joint detection, multiple 
communications are detected simultaneously. 

To utilize these detection techniques, it is desirable to have an estimation of the 
wireless channel in which each communication travels. In a typical TDD system, the channel 
estimation is performed using midamble sequences in communication bursts. 

A typical communication burst 16 has a midamble 20, a guard period 18 and two data 
bursts 22, 24, as shown in Figure 2. The midamble 20 separates the two data bursts 22, 24 
and the guard period 18 separates the communication bursts 16 to allow for the difference 
in arrival times of bursts 16 transmitted from different transmitters. The two data bursts 22, 
24 contain the communication burst's data. The midamble 20 contains a training sequence 
for use in channel estimation. 

After a receiver receives a communication burst 16, it estimates the channel using the 
received midamble sequence. When a receiver receives multiple bursts 16 in a time slot, it 
typically estimates the channel for each burst 16. One approach for such channel estimation 
for communication bursts 16 sent through multiple channels is a Steiner Channel Estimator. 
Steiner Channel Estimation is typically used for uplink communications from multiple UEs, 
14 t to 14 3 , where the channel estimator needs to estimate multiple channels. 

In some situations, multiple bursts 1 6 experience the same wireless channel. One case 
is a high data rate service, such as a 2 megabits per second (Mbps) service. In such a system, 
a transmitter may transmit multiple bursts in a single time slot. Steiner estimation can be 
applied in such a case by averaging the estimated channel responses from all the bursts 16. 
However, this approach has a high complexity. Accordingly, it is desirable to have alternate 
approaches to channel estimation. 

SUMMARY 

A single transmitter transmits K communication bursts in a shared spectrum in a time 
slot of a time division duplex communication system. The system associated with N 
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midamble sequences. Each burst has an associated midamble sequence. A receiver receives 
a vector corresponding to the transmitted midamble sequences of the K communication 
bursts. A matrix having N identical right circulant matrix blocks is constructed based in part 
on the known N midamble sequences. The wireless channel between the transmitter and 
receiver is estimated based on in part one of the N blocks and the received vector. 

BRIEF DESCRIPTION OF THE DRAWING(S) 
Figure 1 is a wireless communication system. 
Figure 2 is an illustration of a communication burst. 
Figure 3 is a simplified multiburst transmitter and receiver. 
Figure 4 is a flow chart of multiburst channel estimation. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 
Figure 3 illustrates a simplified multicode transmitter 26 and receiver 28 in a 
TDD/CDMA communication system. In a preferred application, such as a 2 Mbs downlink 
service, the receiver 28 is in a UE 14 x and the transmitter 26 is in a base station 12 l5 although 
the receiver 28 and transmitter 26 may be used in other applications. 

The transmitter 26 sends data over a wireless radio channel 30. The data is sent in K 
communication bursts. Data generators 32 x to 32 K in the transmitter 26 generate data to be 
communicated to the receiver 28. Modulation/spreading and training sequence insertion 
devices 34! to 34 K spread the data and make the spread reference data time-multiplexed with 
a midamble training sequence in the appropriate assigned time slot and codes for spreading 
the data, producing the K communication bursts. Typical values of K for a base station 12 1 
transmitting downlink bursts are from 1 to 16. The communication bursts are combined by 
a combiner 48 and modulated by a modulator 36 to radio frequency (RF). An antenna 38 
radiates the RF signal through the wireless radio channel 30 to an antenna 40 of the receiver 
28. The type of modulation used for the transmitted communication can be any of those 
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known to those skilled in the art, such as binary phase shift keying (BPSK) or quadrature 
phase shift keying (QPSK). 

The antenna 40 of the receiver 28 receives various radio frequency signals. The 
received signals are demodulated by a demodulator 42 to produce a baseband signal. The 
5 baseband signal is processed, such as by a channel estimation device 44 and a data detection 
device 46, in the time slot and with the appropriate codes assigned to the transmitted 
communication bursts. The data detection device 46 may be a multiuser detector or a single 
user detector. The channel estimation device 44 uses the midamble training sequence 
component in the baseband signal to provide channel information, such as channel impulse 
10 responses. The channel information is used by the data detection device 46 to estimate the 
transmitted data of the received communication bursts as hard symbols. 

O 

:^ To illustrate one implementation of multiburst channel estimation, the following 

§ 1 midamble type is used, although multiburst channel estimation is applicable to other 

5= midamble types. The K midamble codes, , where k = 1 ...K, are derived as time shifted 

S-.J 

QB5 versions of a periodic single basic midamble code, m p , of period P chips. The length of each 
midamble code is l m = P + W -1 . W is the length of the user channel impulse response. 

: s 

tii Typical values for L m are 256 and 512 chips. W is the length of the user channel impulse 

U response. Although the following discussion is based on each burst having a different 
midamble code, some midambles may have the same code. As, a result, the analysis is based 

20 on N midamble codes, N < K. Additionally, the system may have a maximum number of 
acceptable midamble codes N. The receiver 28 in such a system may estimate the channel 
for the N maximum number of codes, even if less than N codes are transmitted. 

The elements of m p take values from the integer set { 1 , - 1 } . The sequence m p is first 

converted to a complex sequence m P [/] = / -m ? [/], where / =!•••? . The m [k) are obtained 
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by picking K sub-sequences of length L m from a IP long sequence formed by 
concatenating two periods of m P . The /'* element of m {k) is related to m P by Equation 1. 

m^ = m P [(K~k)W for l<i<P-(K-k)W 

= in P [i-P+(K-k)w\ for P-(K-k)W <i<P+W-l 

Equation 1 

Thus, the starting point of m [k \ k =l...K shifts to the right by W chips as k increases from 
1 to K. 

The combined received midamble sequences are a superposition of the K 
convolutions. The k h convolution represents the convolution of m {k) with ^(*) . ^(*) is the 

channel response of the k* h user. The preceding data field in the burst corrupts the first 
(W - 1) chips of the received midamble. Hence, for the purpose of channel estimation, only 

the last P of L m chips are used to estimate the channel. 

Multiburst channel estimation will be explained in conjunction with the flow chart of 
Figure 4. To solve for the individual channel responses , Equation 2 is used. 



mp 
m\ 


' * rri(K-\)w+\ 


m(K-\)\v 


rri(K-2)w+\ i 

" ffl(K-2)W+2 1 


i raw 


• mi 
- mi 




(0~ 

h 

(2) 

h 




rw 
rw+\ 


mi 


" 1Tl(K-l)W+3 


rri(K-\)w+2 ' 


" rri(K-2)w+3 ! ' 


* ! mw+2 ' 


* TU3 


X 






rw+2 




: rri(K-\w 


rri{K-\)w~\ 


m(K-2)w ! 


! mw-\ 


mp 




h 







Equation 2 




I-2-149.1US 



r w * * ' t lm ^ ^ e received combined chips of the midamble sequences. The m values are the 
elements of m p . 

Equation 2 may also be rewritten in shorthand as Equation 3. 



2M {k) h {k) =r 



Equation 3 



t-i 



Each M {k) is a KW-by-W matrix, r is the received midamble chip responses. When all the 



bursts travel through the same channel, • • -h [k) can be replaced by h as in Equation 4, 
50. 

K 



£m«> h = 



k=\ 



Equation 4 



G is defined as per Equation 5. 

G =[m (1) ,...,M (;:) ,...,M ( ^ ) ] Equations 

As a result, G is a KW-by-KW matrix. Since G is a right circulant matrix, Equation 4 can 
be rewritten using K identical right circulant matrix blocks B, as per Equation 6, 52. 

B 



k=\ 



B 
B 



= D 



Equation 6 



B is a W-by-W right circulant matrix. The number of B-blocks is K. Using Equation 6, 
Equation 4 can be rewritten as Equation 7. 
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Dh =r 



Equation 7 



Equation 7 describes an over-determined system with dimensions KW-by-W. One approach 
to solve Equation 7 is a least squares solution, 54. The least squares solution of Equation 7 
is given by Equation 8. 



5 h=(D H DY l D H r 



Equation 8 



D is the hermitian of D. 

Applying Equation 6 to Equation 8 results in Equation 9. 

( D » B )- = !(,«,)- 



Equation 9 



The received vector r of dimension KW can be decomposed as per Equation 10. 



l&O r = 



u 



Equation 10 



The dimension of r k is W. Substituting Equations 9 and 10 into Equation 8, the least-squares 
solution for the channel coefficients per Equation 11 results. 



h = (5 "fi) -1 B H 



K k=l J 



Equation 1 1 
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r k represents the average of the segments of r . Since B is a square matrix, Equation 1 1 
becomes Equation 12. 



h = B~ l r k 



Equation 12 



Since B is a right circulant matrix and the inverse of a right circulant matrix is also right 
circulant, the channel estimator can be implemented by a single cyclic correlator of 
dimension 57, or by a discrete fourier transform (DFT) solution. 

A W point DFT method is as follows. Since B is right circulant, Equation 13 can be 

used. 



B=D- l -A c -D w 



D w is the W point DFT matrix as per Equation 14. 



Equation 13 
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Equation 14 



A c is a diagonal matrix whose main diagonal is the DFT of the first column of B, as per 
Equation 15. 



A c =diag(D w (B{.,\))) 



Equation 15 
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W = e jW . Thus, D w is the DFT operator so that D w x represents the W point DFT of the 

vector x . By substituting Equation 13 into Equation 12 and using z)" 1 =-^— , results in 

W 



Equation 16. 



h = 



Equation 16 



is the element-by-element complex conjugate of D w . 
Alternately, an equivalent form that expresses jj~ in terms of instead of A c can 
be derived. A^ is a diagonal matrix whose main diagonal is the DFT of the first row of B 
per Equation 17. 

A R =diag(D w (fl(l,:))) Equation 17 

Since the transpose of B, B T , is also right circulant and that its first column is the first row 
of B, B T can be expressed by Equation 18. 

B T =Dy -A R -D W Equation 18 

Using Equation 18 and that = D W ,A? R =A R and that for any invertible matrix A, 

7 ) 1 = (a _1 J , B can be expressed as per Equation 19. 
B=D W -A R - D^ 1 Equation 19 

Substituting Equation 19 into Equation 12 and that D^ 1 = ~~r results in Equation 20. 
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h = 



1 



d w -a-^-d; 

W J 



Equation 20 



Equations 16 or 20 can be used to solve for h . Since all DFTs are of length W, the 
complexity in solving the equations is dramatically reduced. 

An approach using a single cycle correlator is as follows. Since B _1 is the inverse 
of a right circulant matrix, it can be written as Equation 21 . 
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Equation 21 



Tx Tw 
T 2 Tx 



The first row of the matrix T is equal to the inverse DFT of the main diagonal of A^ 1 . Thus, 
the matrix T is completely determined by A^ 1 . 

The taps of the channel response h are obtained successively by an inner product of 
successive rows of T with the average of W- length segments of the received vector r . The 
successive rows of T are circularly right shifted versions of the previous row. Using registers 
to generate the inner product, the first register holds the averaged segments of r , and the 
second register is a shift register that holds the first row of the matrix T. The second register 
is circularly shifted at a certain clock rate. At each cycle of the clock, a new element of h 
is determined by the inner product of the vectors stored in the two registers. It is 
advantageous to shift the first row of the matrix T rather than the received midambles. As 
a result, no extra storage is required for the midambles. The midambles continue to reside 
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in the received buffer that holds the entire burst. Since the correlator length is only W, a 
significant reduction in complexity of estimating the channel is achieved. 

* * * 
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